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pg . The ~ 250 kg highly radiopure Nal(Tl) DAMA/LIBRA apparatus, running at the 

./ I Gran Sasso National Laboratory (LNGS) of the LN.F.N., is described. 
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1 Introduction 

While running the first generation DAMA/Nal apparatus (~ 100 kg highly 
radiopure Nal(Tl)) [Tll2llMll5]l6]l7ll8MT0llTT^ DAMA proposed to realize a 
ton apparatus in 1996 [13]. Thus, a second generation R&D project for highly 
radiopure Nal(Tl) detectors was funded and carried out over several years 
with leader companies in order to realize, as a new step, a second genera- 
tion experiment with an exposed mass of about 250 kg. The ~ 250 kg highly 
radiopure Nal(Tl) DAMA/LIBRA (Large sodium Iodide Bulk for RAre pro- 
cesses) apparatus has been designed and built mainly to further investigate 
the Dark Matter particle component (s) in the galactic halo by exploiting the 



annual modulation signature; it also aims to improve the investigation of sev- 
eral rare processes. DAMA/Nal completed its data taking in July 2002 and 




Fig. 1. Disinstalling DAMA/Nal and installing DAMA/LIBRA; all procedures in- 
volving detectors, photomultipliers (PMTs), etc. have been carried out in HP (High 
Purity) Nitrogen atmosphere. Improving the experimental site (walls, floor, etc) af- 
ter removal of the DAMA/Nal electronics {top left); dismounting the DAMA/Nal 
detectors in HP Nitrogen atmosphere {top right); assembling a DAMA/LIBRA de- 
tector in HP Nitrogen atmosphere {bottom left); installing the DAMA/LIBRA de- 
tectors in HP Nitrogen atmosphere {center right) and a view of the inner Cu box 
filled by the DAMA/LIBRA detectors in HP Nitrogen atmosphere {bottom right). 

the installation of DAMA/LIBRA started (see Fig. [1]). DAMA/LIBRA has 



begun first operations in March 2003. 

In the present paper the DAMA/LIBRA apparatus and, in particular, the 
features of its new detectors, developed by means of the second generation 
R&D, are described. 



2 The target material 

Highly radiopure Nal(Tl) scintillators have been chosen as target - detec- 
tor material for particle Dark Matter investigations since they offer many 
competitive aspects; among them we remind: 1) well known technology; 2) 
high radiopurity reachable by material selections and protocols, by chemi- 
cal/physical purifications, etc.; 3) large mass feasible; 4) high duty cycle; 5) 
routine calibrations feasible down to keV range in the same conditions as the 
production runs; 6) well controlled operational conditions and monitoring fea- 
sible; 7) absence of microphonic noise; 8) suitable signal/noise discrimination 
near the energy threshold profiting of the relatively high available number of 
photoelectrons/keV and of the well different timing structures of the PMT 
noise pulses (single fast photoelectrons with decay time of order of tens of ns) 
with respect to the Nal(Tl) scintillation pulses (time distribution with time 
decay of order of hundreds of ns); 9) high light response, that is keV threshold 
really reachable; 10) no necessity of re-purification or cooling down/warming 
up procedures (implying high reproducibility, high stability, etc.); 11) possi- 
bility to exploit the granularity of the apparatus, an interesting feature for 
Dark Matter particle investigations and for background recognition; 12) sen- 
sitivity to both high (by Iodine target) and low (by Na target) mass WIMP 
(Weakly Interacting Massive Particles) candidates; 13) high sensitivity to the 
class of WIMP candidates with spin-independent (SI), spin- dependent (SD) 
and mixed (SI&SD) couphngs (see e.g. [51,13]); 14) high sensitivity to several 
other existing scenarios (see e.g. [H.5;6;7;8;9, 10,11, 15]) and to many other pos- 
sible candidates including those producing just electromagnetic radiation in 
the interaction; 15) possibility to effectively investigate the annual modula- 
tion signature in all the needed aspects; 16) "ecologically clean" apparatus, 
thus no safety problems; 17) technique cheaper than every other considered in 
the field; 18) small underground space needed; 19) pulse shape discrimination 
feasible at reasonable level when of interest. 

Moreover, highly radiopure Nal(Tl) scintillators can also offer the possibility 
to achieve significant results on several other rare processes - as already done 
by the DAMA/Nal apparatus - such as e.g. [12]: i) possible violation of Pauli 
exclusion principle in ^^^I and in ^^Na; ii) electron stability; iii) charge non- 
conserving processes; iv) search for solar axions; v) search for exotic matter; vi) 
search for spontaneous transition of nuclei to a superdense state; vii) search 
for spontaneous emission of heavy clusters in ^^^I; viii) search for possible 
nucleon, di-nucleon and tri-nucleon decays into invisible channels; etc. 



3 Layout of the experimental apparatus 



The main parts of the experimental apparatus are: i) the installation; ii) the 
multicomponent passive shield; iii) the 25 highly radiopure Nal(Tl) detectors; 
iv) the glove-box for calibrations; v) the electronic chain and the monitor- 
ing/alarm system; vi) the data acquisition (DAQ) system. 

All the materials constituting the apparatus have been selected for low ra- 
dioactivity by measuring samples with low background Ge detector deep un- 
derground in the Gran Sasso National Laboratory and/or by Mass Spectrom- 
etry (MS) and by Atomic Absorption Spectroscopy (AAS). 



3.1 The installation 



A schematic view of the building housing the DAMA/LIBRA experimental 
apparatus is shown in Fig.[2l Even the materials, the paints, etc., used to build 
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Fig. 2. Schematic view of the building hosting the DAMA/LIBRA apparatus (not to 
scale). About 1 m concrete (made from the Gran Sasso rock material) almost fully 
surrounds (mostly outside the building) the apparatus. In particular, the room A 
hosts the glove box for calibrations; some facilities are in the room B] the sensitive 
parts of the apparatus are in the room C inside the multi-component shield as well 
as the electronics; the computer for DAQ is in the room D. 

the installation, have been selected for low radioactivity as much as possible 
at time of construction. Most of the materials have been underground since 
about 15 years or more. 

At the ground floor, the room D in Fig. [2] houses the computer devoted to 
the data acquisition and the room C contains the sensitive part of the ap- 



paratus, the passive shield, the electronic chain and the acquisition system. 
The floor and the walls of the room C are sealed from external air by a tight 
coverage made by Supronyl (permeability equal to 2 x 10~^^ cm^/s ^\)] this 
is the first level of sealing from the Radon present in trace in the air of the 
underground laboratory. The door at entrance of the room C is air tight and 
- since HP Nitrogen is released in this room - a low level oxygen alarm is 
in operation. The level of the environmental Radon inside the experimental 
room C is continuously monitored by a Radon-meter and recorded with the 
production data; it shows that the environmental Radon in the installation is 
at level of the Radon-meter sensitivity (~ 3 Bq/m'^). 

The DAMA/LIBRA passive shield is supported by a metallic structure mounted 
above a concrete basement (see Fig. [2]). A neoprene layer separates the concrete 
basement and the floor of the laboratory. The space between this basement 
and the metallic structure is filled by paraffin for several tens cm in height. 
This shield is sealed in a plexiglass box maintained in HP Nitrogen atmo- 
sphere in slight overpressure with respect to its external environment; this is 
the second level of sealing from the external air. 

Two mechanical systems allow the lowering of the front and the back sides 
of the shield through an electric engine and, hence, to access the detectors 
during their installation (also see Fig. [1]). 

On the top of the shield (at the first fioor, room A) a glove-box (also con- 
tinuously maintained in the HP Nitrogen atmosphere) is directly connected 
through Cu pipes to the inner Cu box, housing the detectors. The inner Cu 
box is also continuously fiushed with HP Nitrogen and maintained in slight 
overpressure with respect to its external environment (this is the third level 
of sealing from external air in the underground laboratory). The pipes are 
filled with low radioactivity Cu bars (covered by 10 cm of low radioactive Cu 
and 15 cm of low radioactive Pb) which can be removed to allow the inser- 
tion of radioactive sources holders for calibrating the detectors in the same 
running condition without any contact with the installation environment. The 
glove-box is also equipped with a compensation chamber. 

Paraffin/polyethylene fill the space surrounding the plexiglass box as much 
as possible depending on the available space. Moreover, mostly outside the 
building, the DAMA/LIBRA apparatus is almost fully surrounded by about 
1 m concrete made from the Gran Sasso rock material. 

The whole installation is under air conditioning to assure a suitable and stable 
working temperature for the electronics. The huge heat capacity of the multi- 
tons passive shield (~ 10^ cal/°C) assures further a relevant stability of the 
detectors' operating temperature. In particular, two independent systems of 
air conditioning are available for redundancy: one cooled by water refrigerated 
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Fig. 3. A schematic and simplified view of the monitoring/alarm system. Legend: 
1) temperature probe of the detectors' environment; 2) temperature probes of the 
electronic devices; 3) temperature probe of the DAQ system; 4) interface of the 
probes; 5) HP N2 gas inlet in the inner CU box; 6) HP N2 gas inlet in the external 
plexiglass box; 7) and 9) fluxmeters; 8) and 10) valves for the flux regulating; 11) 
interface towards the LNGS alarm system; 12) fan coil; 13) temperature probe for 
the water in the cooling system; 14) cooling system towards a chiller; 15) gas outlet; 
16) pressure gauge; 17) gas bubbler; 18) Radon meter; 19) interface of the control 
of the gas bottles; 20) and 21) pressure meters; 22) and 23) electrovalves; 24) local 
junction box of the safety and alarm system of LNGS. Moreover, self-controlled 
computer processes automatically monitor also several other parameters (the above 
mentioned and those from the DAQ) and manage software alarms (by e-mails and 
by SMS) to operator. 



by a devoted chiller and the other operating with cooling gas. 



A hardware/software monitoring system provides data on the operating condi- 
tions. In particular, several probes are read out and the results are stored with 
the production data. A schematic view of the monitoring/alarm system is re- 
ported in Fig. [3l It is directly interfaced to the safety system of the laboratory, 
for prompt action of the LNGS operator. Moreover, self-controlled computer 
processes automatically monitor several parameters - including those from 
DAQ - and manage alarms (by e-mails and by SMS) to the DAMA operator 
if one of them would go out of the allowed range (whose restrictions depend 
on the required stability level of the running conditions). 



4 The shield from environmental radioactivity 



Figure UK) shows a schematic view of the Cu/Pb/Cd/polyethylene-paraffin 
low background hard shield against environmental radioactivity made of very 
high radiopure materials, which are underground since at least about 15 years. 
Moreover, paraffin/polyethylene fill the space surrounding the plexiglass box 
as much as possible depending on the available space and, as mentioned, 
mostly outside the building, the DAMA/LIBRA apparatus is almost fully 
surrounded by about 1 m concrete, made from the Gran Sasso rock, which 
acts as a further neutron moderatoj ^ I. 

The 25 highly radiopure Nal(Tl) detectors are placed in 5 rows by 5 columns. 
Each one has 9.70 kg mass and the size is: 10.2 x 10.2 x 25.4 cm^. They 
constitute the sensitive part of the DAMA/LIBRA apparatus and together 
with the PMTs (two for each detector at opposite sides) and their Cu shields 
are enclosed in a sealed low radioactive OFHC Cu box continuously flushed 
with HP N2 longly stored deep underground (see Fig. [5]). The Cu box is 
maintained at small overpressure with the respect to the environment, such 
as also the glove-box on the upper level (see Fig. [2]). 

In the DAMA/LIBRA apparatus the PMTs and their 10 cm long light guides 
(made of Suprasil B [22] and directly coupled to the bare crystal acting also 
as optical windows) are surrounded by shaped low-radioactive OFHC freshly 
electrolyzed copper shields; see Fig. Wp). 

All the materials constituting the Cu box and the Cu bricks, used to fill it as 
much as possible, have been selected for low radioactivity. Particular care has 
been devoted to the long-living so-called standard contaminants: ^^'^U, ^^^Th 
and ^°K, since the first two have very rich chains and the latter one is present 
at a level of 0.0117% in the ""*K, which is very abundant in nature. 

The residual radioactivity in some components of the Cu box, that houses the 
Nal(Tl) detectors, is: i) in the copper: < 0.5 ppb, < 1 ppb and < 0.6 ppm of 
238-[j^ 232rj]_^ g^^^ nat-^^ respectively; ii) in the feedthroughs: < 1.6 ppb and < 
1.8 ppm of ^^^Th and ^"-^K, respectively; iii) in the neoprene used in 0-ring: 
< 54 ppb and < 89 ppm of ^'^^Th and "'^*K, respectively. All the limits - at 
the level of sensitivity of the used Ge detector - are at 95% C.L.. 

Outside the Cu box, the passive shield is made by > 10 cm of OFHC low ra- 
dioactive copper, 15 cm of low radioactive lead, 1.5 mm of cadmium and about 
10-40 cm of polyethylene/paraffin (thickness fixed by the available space); see 
Fig. UK). The residual contaminants in some components of the passive shield 



^ Neutron fluxes measured deep underground in the Gran Sasso National Labora- 
tory in the various energy regions are reported in refs. [18|19|20|21j 
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Fig. 4. a) A schematic view of the passive shield of the DAMA/LIBRA apparatus. 
Mostly outside the installation, the DAMA/LIBRA apparatus is also almost fully 
surrounded by about 1 m concrete made of the Gran Sasso rock, b) Up: a scheme 
of the shaped low-radioactivity copper shield for the PMTs; Down: detectors dur- 
ing installation; in the central and right up detectors the new shaped Cu shield 
surrounding light guides and PMTs was not yet applied. See text. 

are (95% C.L.): i) in the copper see before; ii) in the Boliden Lead: < 8 ppb, < 
0.03 ppb and < 0.06 ppm of ^^^U, ^^^Th and ""*K, respectively; ill) in the Boh- 
den2 Lead: < 3.6 ppb, < 0.027 ppb and < 0.06 ppm of ^^su, '^^'^Th and "'**K, 
respectively; iv) in the Polish Lead: < 7.4 ppb, < 0.042 ppb and < 0.03 ppm 
of ^^^U, 232rp]^ ^^^ naij^^ respectively; v) in the polyethylene: < 0.3 ppb, < 0.7 
ppb and < 2 ppm of ^^*U, ^^^Th and ""^^K, respectively; vi) in the plexiglass: 
< 0.64 ppb, < 27.2 ppb and < 3.3 ppm of ^^^U, ^^^Th and ""*K, respectively. 

As mentioned, a sealed plexiglass box (also continuously flushed with HP Ni- 
trogen gas) encloses the passive shield. The Supronyl coverage (see above) and 
the mentioned about 1 m concrete complete the DAMA/LIBRA shield. The 
sealed Cu box, the plexiglass box (both flushed with HP N2) and the Supronyl 
represent the 3-level sealing system from environmental air in the laboratory. 
Moreover, the Radon level inside the installation (that is just after the first 
more external level of sealing from external air) is continuously monitored and 
it is at the level of sensitivity of the Radon-meter; it is continuously recorded 
with the production data. 

Specific procedures for preparing and handling the lead and copper bricks have 
also been selected. In particular, they have been chemically etched in a class 
1000 clean room by using, respectively, HNO3 and HCl aqueous solution in 
iper-pure water (conductivity = 0.0548 /iS/cm at 25 °C). As an example, the 
main steps followed for the inner Cu etching are summarized here: i) Vessel 
I: pre-washing of the brick in iper-pure water; ii) Vessel II: washing in 1.5 1 of 




a) 

Fig. 5. a) Sealing the low radioactive OFHC Cu box. b) The space in front of the 
Cu box is filled by low radioactive Cu bricks and the front wall of the shield can go 
up. 

super-pure HCl 3M; iii) Vessel III: first rinse with iper-pure water (bath); iv) 
Vessel IV: second rinse with iper-pure water (current); v) Vessel V: washing 
in 1.5 1 of ultra-pure HCl 0.5M; vi) Vessel VI: first rinse with iper-pure water 
(bath); vii) Vessel VII: second rinse with iper-pure water (current); viii) Ves- 
sel VIII: third rinse with iper-pure water (current); ix) The Cu brick is dried 
with selected clean towels and HP N2 flux; x) Finally, the Cu brick is sealed 
in two envelopes (one inside the other) previously flushed and filled with HP 
N2 and stored underground until their mounting. The residual contaminants 
in the HCl used in solution with iper-pure water are certified by the producer. 
In particular, standard contaminants are quoted to be: 10 ppb of ""^^K and 
1 ppb of U/Th for super-pure HCl and 100 ppt of """"^K and 1 ppt of U/Th 
for ultra-pure HCl. For each brick the bath was changed and after each step 
the solution of the bath was analyzed with ICP-MS technique. Residual con- 
taminants were preliminarily checked in order to optimize the choice of the 
materials (in particular for gloves) and the cleaning procedure. Moreover spe- 
cial tools have been built (by using teflon and high radiopure OFHC copper) 
to manage the bricks in order to minimize the contact with gloves in presence 
of acid solution. ICP-MS measurements assure the same radiopurity charac- 
teristics of the iper-pure water used in the last rinse, before and after the 
rinse. 



5 The new highly radiopure Nal(Tl) detectors 



The ^ 250 kg Nal(Tl) DAMA/LIBRA apparatus uses 25 Nal(Tl) highly ra- 
diopure detectors with 9.70 kg mass each one (10.2 x 10.2 x 25.4 cm^ volume) 
placed in five rows by five columns. The granularity of the apparatus is an 
interesting feature for Dark Matter particle investigations since Dark Matter 
particles can just contribute to events where only one of the 25 detectors fires 
{single — hit events) and not to whose where more than one detector fire in 
coincidence {multiple — hit events). 



The new DAMA/LIBRA detectors have been built by Saint Gobain Crys- 
tals and Detectors company. The constituting materials have been selected 
by using several techniques; moreover, chemical/physical purifications of the 
selected powders have been exploited. In addition, the more effective growing 
procedure has been applied and new rules for handling the bare crystals have 
been fixed. Each detector is sealed in low radioactivity freshly electrolyzed 
copper housing and has two 10 cm long highly radiopure quartz (Suprasil B) 
light guides which also act as optical windows being directly coupled to the 
bare crystals. 

The bare crystals have been built by Kyropoulos method [23] in a platinum 
crucible according to a strict growing and handling protocol developed and 
agreed under confidential restriction ^ I. The seed, used in the growing, was 
from previous growth for highly radiopure Nal(Tl). The Kyropoulos method 
- following the specific protocol agreed with the Saint Gobain Crystals and De- 
tectors company - has been adopted since it assures in particular: i) possibility 
to operate in a well controlled overall situation for (radio) purity; ii) significant 
(radio)purification during the growing process; iii) high performances of the 
detectors. Other growing processes for Nal(Tl), such as the Bridgeman [23] 
was tested in the past and abandoned since it offered in our experience lower 
performances, mainly as regards the (radio) purification during growth. 

The used Nal powders have preliminarily been selected for radiopurity by var- 
ious methods; then, chemical/physical procedures for the further radiopurifi- 
cation have been exploited by Saint Gobain Crystals and Detectors company. 
The final residual contaminants in the powders are 0.02 ppb for ^^^U, 0.02 
ppb for ^^^Th and < 0.1 ppm (95% C.L.) for "■"■^K, respectively, as measured 
by MS and AAS. As regards the selected Til powder, the residual contami- 
nants in the powders after the chemical/physical purifications are 0.8 ppb for 
238U, 1.2 ppb for 232Th and < 0.06 ppm (95% C.L.) for ""*K, respectively, as 
measured by MS and AAS. It is worth noting that the concentration of Tl in 
the detectors is about 0.1%. Moreover, the growth with Kyropoulos method in 
platinum crucible with a peculiar devoted protocol acts as a further relevant 
radiopurification step. Obviously, the maximal sensitivity to determine the fi- 
nal residual contaminations is reached when directly measured in the crystals 
themselves (see later). In fact, the AAS and MS are mainly limited by the 
radiopurity of the line and by the added materials and - as mentioned - the 
used crystallization process itself acts as an additional considerable purifica- 
tion step. However, it is worth noting that the reachable final levels depend 
both on the initial impurity concentrations and on the used processes and 
handling protocols. Particular care has been devoted to the selection of tools 
and abrasives to be used for cutting and polishing the crystals from the pro- 



^ The confidential agreement protects the intellectual properties both of the com- 
pany and of the collaboration on the specific technical details. 



10 



duced large ingots; this requires, in particular, severe selections and control of 
tools and of procedures in order to avoid any significant pollution on surface 
during the operation. 

The number of materials entering in the production and in the assembling 
of a detector - after the ingots production - has been minimized and deeply 
selected for radiopurity. In particular selected Tetratec-tefion foils, wrapping 
the bare cystals, have been used as diffusers of the scintillation light because 
of their suitable performances and mainly of their radiopurity with respect to 
other possibilities. For the same reason, OFHC freshly electrolyzed copper has 
been used for the housings of the detectors. Preliminary suitable etching of 
the copper housings has been performed for a deep cleaning of possible resid- 
uals from machinery. Obviously, the foam (not distinguishable from outside 
by eyes), usually used inside the housing of standard Nal(Tl) detectors for 
compensation, is not used in the construction of the DAMA/LIBRA detectors 
because of the low background requirements. Instead, a mechanical compen- 
sation for the different specific heats of the Cu and Nal(Tl) is included in the 
design of the copper housing to avoid any possible damages during transporta- 
tion in non-air-conditioned environment before installation deep underground. 
Many other specific expedients have also been exploited. In particular we re- 
mind that the detectors, built for the DAMA apparata in Gran Sasso, have 
always been cut and assembled just after the growth of the crystalline bulk 
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Fig. 6. First moment, r, of the time distribution for each event (recorded by the 
Waveform Analyser) - calculated within a 600 ns averaging time - as a function 
of the energy in unit of keV electron equivalent for four DAMA/LIBRA detectors. 
The two populations 7/e and q's are clearly separated; the a's have shorter r's. 
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in a glove-box in controlled atmosphere, tightly sealed according to special 
requirements, road transported and immediately brought underground. 

A preliminary estimate of the residual U and Th contamination in a detec- 
tor can be obtained from the internal a particles produced by both chains. 
They can be measured by exploiting the a/e pulse shape discrimination in 
Nal(Tl), which has practically 100% effectiveness in the MeV range. Fig. [H] 
shows some examples of the distribution of r as a function of energy for four 
DAMA/LIBRA detectors. The variable r = ^^ is the first moment of the 
time distribution of each scintillation pulse (as recorded by the Waveform 
Analyser; see Sect. [T]); there hi is the pulse height at the time ti and the sum 
is over 600 ns after the starting of the pulse. The two populations 7/e and a's 
are clearly separated; the a particles have shorter r values. The measured a 
yield in the DAMA/LIBRA detectors ranges from 7 to some tens a/kg/day. 

The time-amplitude method'^ I has been used to determine the activity of the 
^^^Th (Ti/2 = 1.912 yr) subchain in the ^^^Th family. In fact, the arrival time 
and energy of the events in each crystal can be used to identify a sequence of 
fast a decays following 228Th decay: ^s^Ra (Q^ = 5.789 MeV, T1/2 = 3.66 d) 
-^ 220r^ ^q^ _ 5405 MeV, T1/2 = 55.6 s) -^ ^lep^ (q^ = 5,905 MeV, T1/2 
= 0.145 s) ^ 212p^_ 

An example of the obtained results for the detector d) of Fig. [H] is shown in 
Fig. [71 The subchain has been selected by looking for events in the energy 
window 3050-3900 keV, followed within 60 s by another event in the energy 
window 3400-4500 keV, still followed within 1 s by the last event in the energy 
window 3650-5100 keV. 

The a peaks as well as the distributions of the time intervals between the 
events are in a good agreement with those expected for the a's from the ^^^Ra 
^ 220Rn ^ 2i6po ^ 2i2pb chain. 

Moreover, the position of the a peaks, reported on an energy scale calibrated 
with 7 sources, allows the determinination of the a//3 light ratio in the used 
crystal; it is: a//3 = 0.467(6) + 0.0257(10) x E^\MeV\, where E^\MeV\ is the 
energy of the a particle in MeV. 

The number of the triple delayed coincidences, selected by the given procedure, 
is 3310 in an exposure of 8100 kg x day. It corresponds, after applying the 



^ The technique of the time-amplitude analysis is based on the identification of 
radioactive isotopes by means of the characteristic time sequence of peculiar events 
(q, /3, 7 of known energy) given by the subsequent decays of daughter nuclides; see 
e.g. ref. 
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Fig. 7. Top: a peaks of ^^^Ra, ^^''Rn and ^^^Po selected by the time-amplitude 
analysis (see text) from the data accumulated during an exposure of 8100 kg x 
day for the triple delayed coincidence with the detector d) of Fig. [6l The energy 
resolutions (a) are: 75, 83, and 90 keV, well compatible with those expected by 
7's calibrations (see Sect. [8]). Bottom: distributions of the time intervals between 
events due to ^^^Ra and ^^'^Rn a decays (left) and between events due to ^^°Rn 
and ^^^Po a decays {right). The obtained half-lives of ^■^''Rn: (67 it 7) s, and of 
^^^Po: (0.142 lb 0.003) s, are in agreement with the table values (55.6 s and 0.145 s, 
respectively [25]). 

efficiencies, to (9.0 ± 0.4) /xBq/kg. 

Repeating the same analysis for all the crystals of the DAMA/LIBRA appara- 
tus we obtain a ^^^Th activity ranging from 2 to about 30 /xBq/kg depending 
on the crystal. 

For completeness we note that the activity value determined in this way for 
each detector is well consistent with that obtained in an independent analysis 
of the so-called Bi-Po events. The Bi-Po events are due to the ^^^Bi j3 decay 
to 2i2po and to the subsequent ^ispo a (Q« = 8.954 MeV) decay to ^"^Pb 
with Ti/2 = 299 ns q|. They can be clearly identified by the recorded pulse 



^ Also the 238^ radioactive chain produces Bi-Po events, due to ^^^Bi f3 decay 
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shapes over a 2048 ns time window. In particular, for the detector d) of Fig. 
|6] the inferred activity for ^^^Th decay subchain is (9.4 ± 1.5) /iBq/kg, well 
compatible with the previous determination. 

If the whole ^^^Th radioactive chain is assumed at equilibrium (as it can be 
the case considering that the longest half-life of the daughter isotopes, ^^^Ra, 
is 5.75 years), the ^^^Th contents in the crystals typically range from 0.5 ppt 
to 7.5 ppt. 



238 -y- 

An estimate of the ^ssu contamination can be obtained for each detector from 
the measured a activity and the '^^'^Th content evaluated above. Assuming 
- as a first approximation - that also the ^ss-jj ^j^ajn is in equilibrium, one 
obtains values ranging from 0.7 to 10 ppt depending on the crystal. 

However, the hypothesis of equilibrium for the ^ss^j Qj^ain in the detectors is 
not confirmed by the study of the energy distributions of the a particles, which 
can allow in principle the determination of the various contributions from the 
^'^'^U subchains. In Fig. [8] we show - as an example - the distributions of a 
from ^^^U and ^^^Th chains in some of the new Nal(Tl) crystals as collected 
in a live time of 570 hours; there the a energies are given in keV electron 
equivalent. 

In particular, starting from the low energy peak, the five a peaks in the energy 
spectra of Fig. [8] can be associated witlQ: i) ^^^Th (Q„ = 4.08 MeV) + ^^^V 
(4.27 MeV); ii) ^s^u (4.86 MeV) + ^^^Th (4.77 MeV) + ^^^Ra (4.87 MeV); 
iii) 2iopo (5.41 MeV) + ^ssTh (5.52 MeV) + 222Rn (5.59 MeV) + 224Ra (5.79 
MeV); iv) ^ispo (6.12 MeV) + 2i2Bi (6.21 MeV) + 220Rn (6.41 MeV) and v) 
2^^Po (6.91 MeV). Thus, the contribution of each a decay has been simulated 
and fitted to the experimental energy spectra (some examples are given in 
Fig. E]) considering the 238u radioactive chain split into five segments (2^^U -^ 
234U ^ 230rjj^ ^ 226j^g^ ^ 2iopb ^ 206pb) ^^^ ^^^ 232rjj^ ^-^^^^ ^^ equilibrium. 

The fit of the measured alpha spectra allows the determination of the activities 
of the five 238-[j gubchains and of the 232'jh chain. The results confirm the 
hypotheses that the 238-[j ^j^ain is broken in these Nal(Tl) crystals. As an 



to 2i4po and to the subsequent 2i4po q, (Q^ = 7.833 MeV) decay to 2i0pb with 
Ti/2 = 164.3/xs. In particular, their rate is connected with the activity of the 226 j^^ 
— > 2i0pb decay chain segment. 

^ It is worth noting that the a associated with the decays of ^^^Po and of 2^^Po are 
not present in the shown a plots because they belong to a Bi-Po event and they are 
mainly vetoed by the acquisition system (see later). 
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Fig. 8. The a energy distributions in some of the Nal(Tl) crystals corresponding to 
a live time of 570 hours. The energy is given in keV electron equivalent. 

example in the detector d) of Fig. [6] and Fig. [8] the ^^^Th and ^ss-q contents 
obtained by the fit are: 1) (8.5 ± 0.5) yuBq/kg of ^^^Th [that is, (2.1 ± 0.1) 
ppt, value in agreement with the two determinations given above using the 
time-amplitude and the Bi-Po analyses]; 2) (4.4±0.7) yuBq/kg for ^^^U — > ^^^U 
decay subchain [that is, (0.35 ± 0.06) ppt of ^^^V]; 3) (15.8 ± 1.6) /iBq/kg for 
234u ^ zsorpj^ _^ 230r|nj^ _^ 226j^g^ ^^^.g^y g^bchaius (they all contribute to the 

same peak); 4) (21.7 ± 1.1) /iBq/kg for ^^^Ra — > ^^°Pb decay subchain and 5) 
(24.2 ± 1.6) /iBq/kg for ^lopb -^ ^oep^ ^jecay subchain. 

As it is clear e.g. from Fig. [HI the residual contaminants may be slightly dif- 
ferent even among detectors made from Nal(Tl) crystals grown with the same 
selection of materials, purification processes and protocols. In fact, some casual 
pollutions during the growth and handling procedures may in principle be pos- 
sible, being the detectors built in an industrial environment. Differences may 
also arise depending on the use of different bulks or on which part of a crystal- 
lized bulk has been used to build the detector. In fact, the purification during 
crystallization may be not uniform in the whole bulk mass. Moreover, the uni- 
formity of the contaminants distribution inside the total material needed to 
construct each part of the detectors cannot be assured. Obviously, casual pol- 
lution may also occur when handling the detectors in industrial environment 
or deep underground without the needed extreme care. 



nat 



K 



An estimate of the potassium content in the DAMA/LIBRA crystals has been 
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obtained investigating over large exposure the presence of peculiar double 
coincidences. In fact, the ^°K (0.0117% of '^'^^K) also decays by EC to the 1461 
keV level of ^°Ar (b.r. 10.66%) followed by X-rays/Auger electrons, that are 
contained in the crystal with efficiency ~ 1, and a 1461 keV de-excitation 7. 
The latter one can escape from one detector (hereafter A) and hit another 
one, causing the double coincidence. The X-rays/Auger electrons give rise in 
the detector A to a 3.2 keV peak, binding energy of shell K in ^"AiL^. 



I 




taw 1400 



Fig. 9. Example of the analysis to determine "-"^K contamination in one of the 25 
crystals (see text). The scatter plot shows the low energy region of the considered 
crystal, A, as a function of the energy detected in the other crystal involved in the 
double coincidence. The threshold of each PMT is at single photoelectron level. For 
comparison, the software energy threshold used in the data analyses of the single-hit 
events for Dark Matter particle investigation: 2 keV, is shown as continuous line. 

The experimental data have been analyzed searching for these double coin- 
cidences; Fig. [9] shows as an example a scatter plot of the energies of the 
detector A and all the other detectors involved in the coincidence. It is evi- 
dent a spot that correlates the 1461 keV events in the other crystals with the 
3.2 keV peak in crystal A. The detection efficiency for such coincidences has 
been evaluated for each crystal by MonteCarlo code. The analysis has given 
for the "■"*K content in the crystals values not exceeding about 20 ppb. It is 
worth noting that the identification of the 3.2 keV peak offers also an intrinsic 
calibration of each DAMA/LIBRA detector in the lowest energy region (see 
Sect. ED. 

125j 

As regards possible non-standard contaminants, we remind that during the 
period of powders' storage and/or crystal growing and handling at sea level 
(of order of at least some months for large crystals) radioactive isotopes can 
be produced by cosmic-ray interactions. Therefore, it is good practice in this 
kind of experiments to wait for the decay of the short-life cosmogenic isotopes 
before considering the production data in underground. 



^ In the (76.3 it 0.2)% of the cases an electron from shell K (E/^ = 3.2 keV) is 
involved in the process, in the (20.9 it 0.1)% an electron from shell L (E^ = 0.3 
keV) and in (2.74 it 0.02)% electron from upper shells. 
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In particular, here we comment the case of cosmogenic ^^^I which decays by 
EC to ^^^Te giving 35.5 keV 7's (or internal conversion electrons) and, at the 
same time, tellurium X-rays and Auger electrons, whose total energy is 31.8 
keV, giving rise to a peak around ~ 67 keV with half-life 59.40 days (see 
Fig. [TOl) . In this case the best situation is reached after ~ 8 months of storage 
deep underground. 
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Fig. 10. a) Difference between the energy distributions collected at Gran Sasso with 
a DAMA/LIBRA Nal(Tl) detector just stored underground and 15 months later. 
The role of ~ 67 keV peak due to ^^^I activation (half-life 59.40 days) is evident. 
In particular, the EC of ^■^^I from shell K (80.12%) gives rise to the 67.3 keV peak, 
while EC from upper shells produces the ~ 40 keV peak, b) The same as a) but 
obtained for a detector stored underground for a longer time; the peaks due to ^^^I 
are not present. 



Some arguments on cosmogenic activation in Nal(Tl) can be found in ref. 
while the production rate of cosmogenic isotopes can be calculated mainly ac- 
cording to the model of ref. |27j. However, the activity of cosmogenic isotopes 
in Nal(Tl) detectors would be in practice generally lower than the calculated 
one because of purification and saturation processes during the growing pro- 
cedures. 



129 



I and ^i"Pb 



The cosmogenic ^^^I (T1/2 = 1.57 x 10^ yr) can be present in the natural Iodine 
with a percentage of the order of 1.5 x 10^^^ (considering sediments below the 
zone of bioturbation, i.e. without the presence of anthropogenic Iodine [28]). 
It decays by /?~ into the excited level of ^^^Xe at 39.57 keV. Thus, this pro- 
cess produces a de-excitation 7 of 39.57 keV (or internal conversion electron) 
plus an electron with energy distribution according to the beta spectrum of 
129 J (ixiaximum energy 154 keV). Therefore, the total energy deposited in the 
crystal from this process has a peculiar distribution with sharp rise around ~ 
40 keV and smooth behaviour up to ~ 194 keV (sum of the gamma energy and 
of the maximum energy of the electrons). Fig. [TT] shows the energy spectra of 
two detectors in the energy region of interest for this process and the model 
of background. Hence, the amount of cosmogenic ^^^I has been estimated to 
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be there i29j/natj ^ (i.7±o.l) x IQ-^^. If this value is used for dating the NaT 
powders, one obtains that they have been extracted from an ore with an age 
of order of 50 Myr. 
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Fig. 11. Energy distributions of two new DAMA/LIBRA Nal(Tl) detectors in 
the energy region of interest for the cosniogenic ^^^I decay process (histogram) 
and the model of background (Hues). The contribution of internal ^'^'^Pb is visible 
only in the plot of the second crystal (< 12 //Bq/kg at 90% C.L. and (29 it 3) 
/iBq/kg, respectively). The amount of cosmogenic ^^^I has been estimated to be 
i29j/natj ^ (^^^ 7 _^ qjJ ^ ^^q-is f^^. ^^^^ ^Yid crystals. 



Fig. [TT] also allows an independent evaluation of the content of ^^°Pb in the 
crystals. The ^^°Pb isotope (3 decays with a branching ratio of 84% into the 
excited level of ^^°Bi at 46.5 keV; the electron has an energy distribution 
according to the beta spectrum (maximum energy: 16.6 keV). The presence of 
a contamination of ^^''Pb mainly on the housings' surfaces has been used by the 
DAMA/Nal apparatus as a calibration point at 46.5 keV (due to the gamma 
from the de-excitation of the ^^°Bi excited level) during its data taking (see for 
example [2P] ). Such a peak is not so appreciable in the new DAMA/LIBRA 
detectors; anyhow, the contents of internal ^^°Pb can be determined for these 
detectors by investigating the energy region around ~ 50 keV. For the two 
detectors of Fig. [TT] one gets: (7±4) /xBq/kg (safely corresponding to an upper 
limit of < 12 /^Bq/kg at 90% C.L.) and (29 ± 3) /iBq/kg, respectively. Let us 
note that for the second detector in Fig. [TT] - corresponding to the detector 
d) of Fig. [6] and Fig. [8] - an independent determination of ^^°Pb content has 
been mentioned above when studying the a spectrum in the crystals: the two 
values are well compatible. 



Repeating the same analysis for all the new DAMA/LIBRA detectors, we 
obtain that the amount of cosmogenic ^^^I is at the same level for all the new 
DAMA/LIBRA detectors; in addition, the content of ^^'^Pb in the crystals 
typically ranges: (5 — 30) /iBq/kg. 
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Na 



Among the possible isotopes produced by cosmogenic activation at sea level in 
Nal(Tl) we also investigate the ^^Na (T1/2 = 2.6 yr). The calculated maximum 
rate level is ^ 100 cpd/kg at sea level. An estimate of the ^^Na activity in 
the detectors can be obtained by searching for triple coincidences induced by 
the /3+ decay of '^^^a followed by 1274.6 keV de-excitation 7 (b.r. 90.33%). 
In particular, we have searched for events where the positron and one of the 
two 511 keV annihilation 7 release all their energy in one detector while the 
second 511 keV annihilation 7 and the 1274.6 keV 7 hit other two detectors. 
For example, choosing an energy window of 650-1000 keV for the first detector 
and ±2(7 around 511 keV for the second one, we look for the typical peak at 
1274.6 keV in the third detector. Accounting for the detection efficiency of the 
process, the obtained activity in the detectors ranges from an upper limit of 
15 /iBq/kg (90% CL.) to some tens /xBq/kg depending on their arrival time 
underground. 
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Na 



Another source of internal background is the possible presence of Na. 

In fact, environmental neutrons would induce the reaction ^^Na(n,7)^'^Na with 
0.1 barn cross-section and the reaction ^^Na(n,7)^^'"Na with 0.43 barn cross- 
section [29j. The ^''Na isotope is a /3-emitter (end point equal to 1.391 MeV) 
with two prompt associated 7's (2.754 and 1.369 MeV). On the other hand, 
the ^^'"Na isotope decays 100% of the times in ^^Na by internal transition 
with a 7 of 0.472 MeV. Thus, the presence of ^"^Na has been investigated with 
high sensitivity by looking for triple coincidences induced by a /3 (El) in one 
detector and the two 7's in two adjacent ones (E2, E3). In particular, we con- 
sider: i) 0.57 MeV < El < 1.3 MeV, ii) E2 and E3 inside ±lo" energy windows 
around the photopeak positions. Only one event satisfying the requirements 
has been found in the considered data set (106376 kgxday); it can also be as- 
cribed to side processes (for example to the decay of ^os^^ jj-^ ^j-^g 232rj]-^ chain). 
Safely we consider that the number of detected events is lower than 3.89 (90% 
CL.) following the procedure given in |30j; thus, accounting for the detection 
efficiency of the process, an upper limit of 2409 ^^Na decays (90% CL.) is 
obtained for the considered exposure. It corresponds to an activity < 0.26 
/iBq/kg (90% CL.). Assuming the production of these short-life isotopes as 
fully due to a steady thermal neutron flux, an upper limit on the thermal neu- 
tron flux surviving the multicomponent DAMA/LIBRA shield can be derived 
as: < 1.2 X 10"^ cm'^s'^ (90% CL.). 
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Miscellanea 

Finally, as regards the ^H (T1/2 = 12.3 yrjjj, the calculated maximum rate 
is ~ 20 cpd/kg at sea level. However, as mentioned above, in real conditions 
the ^H activity in a Nal(Tl) detector would be essentially lower than the 
calculated maximum level, because it would be extracted in the process of 
Nal(Tl) purification and growth. In our crystals the ^H content has been 
estimated to be: < 0.9 x 10"^ Bq/kg at 95% C.L.. Moreover, for the ^^Rb 
(Ti/2 = 4.8 X 10^° yr), and for the ^^Kr (T1/2 = 10.7 yr) isotopes we just 
mention the obtained upper limits (95% C.L.): < 3 x 10^^ Bq/kg and < 10~^ 
Bq/kg, respectively. 



6 The photomultipliers 



The photomultipliers used in this experiment have been built by Electron 
Tubes Limited and longly stored underground. The PMTs are made of ultra- 
low background glass. The materials entering in the construction of the PMT 
have been selected by various kind of measurements and, in particular, mea- 
surements on samples have also been carried out with low background DAMA/Ge 
detector deep underground. 




Fig. 12. Left: a DAMA/LIBRA PMT with flying leads directly connected to the 
voltage divider. Right: a voltage divider assembled with miniaturized SMD resistors 
and capacitors mounted on thin teflon socket. 

The PMTs have fiying leads (instead of a hard socket) and are directly con- 
nected to suitable voltage dividers (see Fig. [T2]). The voltage dividers have 



'^ The ^H can be produced by spallation processes induced by cosmic rays on ^^Na 
and ^^'^I. In particular, we recall: n + ^^Na ^'^11 + ^-'^Ne(stable) and n + ^'^'^l — > 



^H + ^^^Te (stable), for neutron with energy > 15, 20 MeV. 
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been assembled deep underground with miniaturized SMD resistors and ca- 
pacitors (also selected for radiopurity by low background Ge detector) and 
mounted on thin teflon sockets (see Fig. [T^ . All the solders were performed 
by using low radioactive Boliden lead and low radioactive resin. 

The voltage dividers have been optimized to reach the best signal/noise ratio; 
a scheme is shown in Fig. 



The applied High Voltage is positive; in this way the cathode is mantained to 
ground in order to avoid any voltage difference with the PMT window, light 
guide and the copper housing. 




Fig. 13. The voltage divider scheme. 

The PMTs work near the maximum high voltage - within the linearity plateau 
- in order to assure a suitable energy threshold and to improve the PMT noise 
rejection near it. 

The PMTs have nine high gain, high stability dynodes of linear focused design, 
high quantum efficiency (~ 30% at 380 nm), good pulse height resolution for 
single photoelectron pulses (peak/ valley > 2), low dark noise rate (^ 0.1 kHz), 
and a gain of ~ 10^. In the experiment each detector has two PMTs both to 
increase the light collection and to reduce the PMT noise contribution near 
energy threshold. They work in coincidence at single photoelectron threshold. 

The residual contaminants in the PMTs have been measured with and without 
a bent light guide (i.e. shielding and not shielding the PMT contribution) and 
a low background Nal(Tl); typically ^'^^Th and ^ss-q ^j.^ Jower than few tens 
ppb and "'**K is of order of few tens ppm. The PMT contribution to the 
counting rate is further reduced by using 3" diameter 10 cm long UV light 
guides (made of Suprasil B) which also act as optical windows, as mentioned 
above. The highly radiopure light guides assure both suitable low radioactivity 
(< 1 ppb of ^^^U, < 2 ppb of ^^^Th, and < 1.3 ppm of ""*K as measured by 
Ge detectors) and very suitable optical characteristics for the Nal(Tl) light 
emission, being transparent down to 300 nm. 

Finally, to further reduce the PMTs contribution to the background, in the 
DAMA/LIBRA apparatus each PMT and related light guide are surrounded 
by a shaped low-radioactivity copper shield (see section H] and Fig. Hb)). 
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7 The electronic chain 



Each crystal of the DAMA/LIBRA apparatus is viewed by two PMTs with 
grounded cathode and supplied by positive high voltage. The HV power supply 
for the PMTs is given by a CAEN multichannel voltage supply with voltage 
stability of 0.1%. 



Figure [H] shows the scheme of the electronic devices of a single detector. In 
particular, the part of the chain regarding the analog signals and their process- 
ing are shown. It is also reported the logic devices providing the trigger of the 
single detector. The signals from each PMT are amplified by a preamplifier 
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Fig. 14. The conceptual scheme of the electronic devices of a single detector. Devices 
processing analog signals and those providing the trigger of the single detector are 
shown. 

having 0-250 MHz bandwidth, a factor 10 gain and a voltage integral linearity 
±0.2%. 

The signal from PMT 2 (see Fig. UM is divided in two branches by a resistive 
passive splitter: 19/20 of the signal is sent to the preamplifier, while the re- 
maining 1/20 - suitably delayed - feeds the CAMAC charge ADC 2 {QADC 
high energy) unit. This latter allows the process of the high energy pulses, 
which instead saturate the remaining part of the electronics. The preamplified 
signals of PMT 1 and PMT 2 are the inputs of linear Fan-in/Fan-out devices 
Fl and F2. They produce three copies of each signal: i) the first copies of 
PMT 1 and PMT 2 signals are sent to two separate channels of a Waveform 
Analyser, which records the signals in a time window of 2048 ns. This is ac- 
complished by using fast VXI Tektronix four-channel TVS641A digitizers with 
a sampling frequency of 1 GSample/s and 250 MHz bandwidth; ii) the second 
copies of PMT 1 and PMT 2 signals are added and sent to the input - after 
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a delay line - of the CAMAC charge ADC 1 {QADC low energy) and to a 
Spectroscopy Amplifier (#i, see later); iii) the last copies are used to generate 
the trigger of the crystal. These latter ones are shaped by Timing Filter Am- 
plifiers {TFAl and TFA2) with an integration time of 50 ns; their outputs 
are discriminated by DISCRl and DISCR2 with the thresholds at the single 
photoelectron level. The coincidence {AND 1) between the two logical NIM 
outputs provides the trigger of the considered detector. 

In order to reject afterglows, Cherenkov pulses in the light guides and Bi-Po 
events, a SOOyUS veto occurs after each event selected by the coincidence ( Gate 
1 in Fig. [13]); this veto introduces a negligible systematic error on the measured 
rate (< 10"^). 

A CAMAC Output Register {OR) provides the signal able to activate the 
coincidence module of each detector; this feature can be used, for example, 
during the calibration runs to enable the detector to be calibrated. 

The outputs of the coincidence device {AND 1) provide: i) the signal for a 
CAMAC scaler to count the number of triggers of each detector. The ratio 
of this number to the elapsed time gives the hardware rate, Rhj of the j— th 
detector; ii) the line used in the main trigger electronic chain (see later and 
Fig. dSj); iii) the line used in the logic system managing the individual triggers 
of the Waveform Analysers (see later and Fig. [T6l) : iv) the line giving the start 
to the Gate 1 which - in addition to the veto of the coincidence AND 1 - feeds 
a CAMAC I/R Pattern Recognition {IR)- The IR allows the identification of 
the detector/detectors involved in the main trigger. 

The electronic devices which generate the main trigger of the acquisition are 
shown in Fig. [T3 The main trigger is provided by the logic OR {OR 1) of 
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Fig. 15. The conceptual scheme of the electronic chain providing the main trigger 
of DAQ. 

all the crystals. The main trigger NIM logic pulses are counted by a Scaler. 
Devoted devices manage the trigger only when the acquisition is ready (see 
Fig. [15]); they issue the trigger of the DAQ. The dead time of the acquisition 
is properly evaluated by using the information from the Scaler. The live time 
is estimated straightforward. 
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When a DAQ trigger occurs, the following logic signals are issued to: i) the 
Gate Generator producing the 600 ns gates for the charge ADCs (QADC); ii) 
the Delay Gate Generator giving the strobe signal to the CAMAC I/R Pattern 
Recognition. This device also generates the LAM in the CAMAC system, and, 
therefore, the interrupt to the CPU of the acquisition computer; iii) the logic 
system managing the individual triggers of the Waveform Analysers. 

The logic system of the individual triggers of the Waveform Analysers has 
been introduced in order to speed up the DAQ. In particular, the trigger 
of the Waveform Analysers (see Fig. [TBI) are issued only if the total energy 
deposited in the detectors is in an energy window suitably chosen. For this 
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Fig. 16. The conceptual scheme of the logic system managing the individual triggers 
of the Waveform Analysers. 

purpose, each line feeds a Spectroscopy Amplifier (5*^4) whose gain is equal- 
ized. The analog sum of the SA signals of all the detectors, whose stability is 
monitored by a peak ADC (not shown in the scheme), feeds a Single Channel 
Analyser (SCA) which allows the selection just of events in the chosen en- 
ergy window. A devoted electronic circuit issues the trigger of the Waveform 
Analyser^^ allocating the detectors producing the main trigger. Summariz- 
ing, an individual trigger to a given Waveform Analyser is issued when the 
following three conditions are fulfilled: i) at least one of its corresponding lines 
has fired; ii) the DAQ trigger is present; iii) the total energy of the events is 
in the chosen energy window. Let us remind that for the events with energy 
outside the chosen energy window (e.g. high energy events) the QADC values 
are acquired anyhow. 

The data acquisition system is made of a Workstation by Compaq with Linux 
Operating System, interfaced with the hardware system through MXI-2 and 
GPIB buses (see Fig. [T7|) . The GPIB bus allows the communication with the 
CAMAC crate housing the QADCs, the scalers and the I/O registers, while 
the MXI-2 bus allows to communicate with the three VXI mainframes, where 
the Waveform Analysers are installed. 



^ We remind that each Waveform Analyser has four channels and, therefore, allo- 
cates the signals of two detectors. 
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Fig. 17. Data acquisition system of DAMA/LIBRA. 

The information collected by DAQ, when a DAQ trigger occurs, are: i) the 
time of occurrence of the event; ii) the pattern of the detectors firing; iii) for 
each fired detector the QADC channels and the pulse profiles of the two PMT's 
of the detector as recorded by the Waveform Analyser, when the conditions 
given above are fulfilled. In particular, the areas of the pulse profiles ( TD 
channel) are proportional to the energy released in the detector (see Sect. |8]). 
Other variables can be constructed by the pulse profiles and they can be used 
for noise rejection (see Sect. [8]). For high energy events the energy can be 
estimated in terms of QADC channels (see Sect. [8]). 

Moreover, let us remind that the DAQ also records - together with the produc- 
tion data - the information coming from the monitoring system, as mentioned 
above (also see Fig. [3]). 

In near future the DAMA/LIBRA DAQ will be improved and new Waveform 
Analysers will be introduced. In particular, in the new DAQ a HP Worksta- 
tion with Intel processor (3.4 GHz) and Red-Hat Linux Operating System will 
be interfaced with the hardware system through MXI-3 and GPIB buses. The 
new MXI-3 bus will communicate with the new Waveform Analysers through a 
CompactPCI mainframe. The new Waveform Analysers will be Quad-Channel 
CompactPCI Digitizer Acqiris DC270 with 250 MHz bandwidth and a sam- 
pling frequency of 1 GSample/s simultaneously on all 4 channels. The data will 
be transferred on fiber optics from the Waveform Analysers to PC (electrically 
decoupled) . 



8 General features 



Uniformity of the light collection 

The absence of dead spaces in the light collection has been carefully investi- 
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gated by performing suitable calibrations in different positions of the detector; 
in fact, by irradiating the whole detector with high-energy 7 sources (e.g., 
^^^Cs) from different positions, no significant variations of the peak position 
and energy resolution have been observed. Moreover, the a peaks at high en- 
ergy (see Fig. [7]) and their energy resolutions are well compatible with those 
expected for 7 calibration (see above) and also support the uniformity of the 
light collection within 0.59{tfJ. 

Photoelectrons/keV 

The number of photoelectrons/keV has been derived for each detector from the 
information collected by the Waveform Analyser over its large time window 
(2048 ns). 



In fact, a clean sample of photoelectrons can be extracted from the end part 
of this time window; there the scintillation pulses are completely ended, while 
afterglow single photoelectron signals can be present. As an example. Fig. [T51 
shows a typical experimental distribution of the area of the single photoelec- 
tron's pulses for a DAMA/LIBRA detector. The relative peak value can be 
compared with the peak position of the distribution of the areas of the pulses 
corresponding to a full energy deposition from the 59.5 keV of the ^^^Am 
source. This procedure permits to obtain the number of photoelectrons/keV 
searched for; they ranges from 5.5 to 7.5 photoelectrons/keV depending on 
the detector. 
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Fig. 18. Typical experimental distribution of the area of the single photoelectron's 
pulses for a DAMA/LIBRA detector, as recorded by the Waveform Analyser over 
2048 ns time window. From this picture 6.7 photoelectrons/keV have been inferred 
for this detector. 



^ For example, from Fig. [7] a) the energy resolution is o" = (75 it 3) keV, while the 
expected one from the fit of Fig. [22] is 72 keV. Hence, with the help of a simulation 
of the light collection, an upper limit of about 0.5% can be set on the disuniformity 
of the light collection. 
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Calibrations at low and high energy 

The low (up to ~ 100 keV) and the high (above ~ 100 keV) energy regions 
have been studied by using various external gamma sources and internal X- 
rays or gamma's. In particular, as shown in Fig. [191 gammas from ^^^Am 
and from ^^^Ba external sources provide 30.4 keV (composite), 59.5 keV and 
81.0 keV peaks; moreover, other calibration points can be obtained by: i) the 
internal X-rays at 3.2 keV (see Fig. [9] and Fig. [T9l) : ii) the gammas and/or X- 
rays and/or Auger electrons due to internal ^^^I (see Fig. [10] a) and Fig. [T9|) . 
giving two structures in the energy spectrum at 40.4 keV (composite peak) 
and at 67.3 keV; iii) the gammas at 39.6 keV plus the /? spectrum due to 
cosmogenic ^^^1 (see Fig. [TTl) . 

In running conditions routine calibrations with ^^^Am sources are regularly 
carried out. 

The linearity of the pulse area as a function of the energy of the calibration 
points and the energy resolution in the low energy region are depicted in 
Fig. [20l The energy resolution behavior for the low energy scale has been 
fitted by a straight line: aiE/E = ,^^^ + fii^E\ the best fit values are 

(Xle = (0.448 ± 0.035) and Ple = (9.1 ± 5.1) x 10"^ 

We remind here that, although the optimization is made for the lowest energy 
region, the data are taken on the full energy scale up to the MeV region by 
means QADC's. Thus, the high energy scale of the QADC has been calibrated 
mainly by using external sources of 7 rays, such as: ^^^Cs (662 keV), ^°Co (1173 
keV and 1332 keV and their sum: 2505 keV) and ^^^Ba (mainly 81 keV and 
356 keV); examples are reported in Fig. [2TJ Moreover, 7 rays of 1461 keV due 
to ^'^K decays, tagged by the 3.2 keV X-rays in an adjacent detector (see also 
Fig. [9]), have been used. 

The linearity of the QADC as a function of the energy and the energy res- 
olution in the high energy region are reported in Fig. [251 The energy reso- 
lution behavior for the high energy scale has been fitted by a straight line: 
aHE/E = ,"^^ -|- Phe] the best fit values are ans = (1.12 ± 0.06) and 

^^E{keV) 

Phe = (17±23) X 10""^. The value ans is larger than the corresponding value 
at low energy since - as discussed previously - the signals for high energy 
events are taken only from one PMT. 

Another method to study the linearity is to use events of double coincidences 
when two detectors are irradiated by 7's of known energy, E-y. As an example. 
Fig. ^3^ left shows the scatter plot of the energies, Ea and Eb, of the two 
detectors A and B when using an external ^^^Cs source placed in between. 
The solid line represents the case Ea + Eb = E^; the data points around this 
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Fig. 19. Examples of energy (pulse area) distributions for: i) X-rays/Auger electrons 
(3.2 keV) due to internal ^''K (top left); ii) gammas and/or X-rays and/or Auger 
electrons - mainly two structures in the energy spectrum at 40.4 keV (composite 
peak) and at 67.3 keV - due to internal ^^^I {top right); iii) gammas from ^^^Am 
external source {bottom left); iv) gammas of 81 keV from ^^^Ba external source 
{bottom right). The lines superimposed to the experimental data have been obtained 
by MonteCarlo simulations, obviously taking into account the geometry and all the 
involved materials. 

line correspond to events where the 662 keV 7 has a Compton back-scattering 
in one detector and the scattered 7 is completely absorbed in the other one. 
Fixing a slice - for example at a fixed Eb value - it is possible to extract 
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Fig. 20. Linearity of the pulse area as a function of the energy and energy resolution 
{cr/E) as a function of E~^''^ for low energy calibrations (internal and external 
radiation) . 
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Fig. 21. Energy distributions for the external sources: ""^^^Cs, 662 keV 7 rays {top 
left); ^'^Co, 1173 keV and 1332 keV 7 rays and their sum: 2505 keV {top right); 
^■^■^Ba, mainly 81 keV and 356 keV 7 rays {bottom left). Energy distribution for 7 
rays of 1461 keV {bottom right) due to ^'^K decays, tagged by the 3.2 keV X-rays in 
an adjacent detector (see also Fig. [9|). The line superimposed to the experimental 
data of the ^^^Ba spectrum has been obtained by MonteCarlo simulations. Here in 
abscissa the QADC channel is reported. 
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Fig. 22. Linearity of the QADC as a function of the energy and energy resolution 
{cr/E) as a function of E^^'"^ in the high energy region. We remind that the signals 
- unlike low energy events - for high energy events are taken only from one PMT. 

the peak position on the Ej^ variable. Applying this procedure to the high 
energy calibration data from various sources other points can be added in the 
linearity plot, as done in Fig. ^^ right. 
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Fig. 23. Left: Example of a scatter plot of the energies, Ej^ and Eb, of two detectors 
A and B when using an external ^^"^Cs source placed in between. The solid line 
represents the case Ea + Eb = 662 keV, energy of the 7 rays. The data points 
around this line correspond to events where the 7 has a Compton back-scattering 
in one detector and the scattered 7 is completely absorbed in the other one. Right: 
Linearity of the QADC as a function of the energy: data points from Fig. [22] (filled 
points) and from tagged 7's by double coincidences as described in the text (squared 
points) . 

Procedures for noise rejection near the energy threshold 

The only data treatment which is performed on the raw data is to eliminate 
obvious noise events (whose number sharply decreases when increasing the 
number of available photoelectrons) near the energy threshold. In particular, 
the DAMA/LIBRA detectors are seen by two PMTs working in coincidence 
and this already strongly reduces the noise near the energy threshold for the 
single-hit events (of interest for Dark Matter particles detection), while ob- 
viously noise is practically absent in multiple-hit events (see e.g. the double 
coincidence case of Fig. [9]), since the probability to have random coincidences 
is negligible (< 3 x lO"*^). 

Moreover, the Nal(Tl) is a suitable detector for an effective noise rejection, 
having no contribution from microphonic noise (as it is instead the case e.g. 
of ionizing and bolometer detectors) and being the physical pulses well dis- 
tinguishable from the noise pulses, essentially PMT noise. This effectiveness 
is higher when a high number of photoelectrons/keV is available and when 
PMTs with suitable performances have been built, as in the present case. 

The noise rejection in the single-hit events is performed by considering that 
the scintillation pulses in Nal(Tl) have a time distribution with decay time of 
order of hundreds of ns (typically ~ 240 ns) , while the noise pulses are single 
fast photoelectrons pulses with decay time of order of tens of ns (see Fig. [2^ . 
The different time characteristics of these signals can be investigated building 
several variables from the information recorded for each pulse by the Wave- 
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Fig. 24. Typical pulse profiles of PMT noise (top) and of scintillation event (bottom) 
with the same area, just above the energy threshold of 2 keV. 



form Analyser over the 2048 ns time window. The variables follow different 
distributions for scintillation and noise pulses. In particular, the fractions of 
the pulse areas evaluated over different time intervals, such as: 



Xi 



Area (from 100 ns to 600 ns) 
Area (from ns to 600 ns) 



Xo 



Area (from ns to 50 ns) 
Area (from ns to 600 ns) 



(1) 



offer a very effective procedure I. In fact, the first variable is distributed 



around for noise events and around 0.7 for scintillation pulses; the second 
variable indeed is distributed around 1 for noise and around 0.25 for scintilla- 
tion event a ^^ I. The values of the two variables quoted above, obtained for each 
event of the production data, can be plotted one vs the other, as in Fig. \25\-left 
panels; the signature of the residual noise above the software energy threshold 
with respect to the scintillation events is evident. 



Thus, suitable acceptance windows can be applied to obtain just scintillation 
events. Figure [25] clearly points out that stringent acceptance windows can 
be used to select a data sample clean from noise; the effectiveness of this 
procedure has already been demonstrated by DAMA/Nal [US]. 



^^ Note that the 600 ns time interval corresponds to the usual hardware gate for 
Nal(Tl) pulses and it is equal to ~ 2.5 times the decay time of a scintillation pulse, 
that is the area (from ns to 600 ns) is proportional to the energy of the event. 
^^ Qualitatively, hypothesizing only one decay constant the expected value for scin- 
tillation events should be 0.64 for the first variable and ~ 0.2 for the second one. 
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Fig. 25. Bidimensional plots in the Xi vs X2 plane. The two plots on the left have 
been obtained for the single-hit production data, which contain both the PMT noise 
and the scintillation events. The two plots on the right have been obtained for events 
from a 7 source of suitable strength; thus just scintillation events are present. The 
plots on the top refer to the 2-4 keV energy interval and those on the bottom to 4-6 
keV. The noise is clearly separated from the scintillation events; see text. 

The overall efficiency for Dark Matter investigation 

As mentioned above, two PMTs on each detector work in coincidence at single 
photoelectron threshold. This - considering also the relatively high available 
number of photoelectrons/keV - assures that the required coincidence of the 
two PMTs does not introduce a very significant hardware cut of the events 
near the energy threshold. In fact, it is ~ 0.8 in the 2-3 keV bin and 1 above 
5 keV for a coincidence time window of 50 ns and practically always equal to 
unity in case of a 100 ns time window. 

The efficiencies related to the used acceptance windows for the single-hit events 
in the lowest energy bins can be determined and properly used in order to 
account for scintillation events, that might be lost by applying the acceptance 
windows. These efficiencies, for each considered energy bin, are measured by 
applying the same acceptance windows to events induced by a ^^^Am source 
of suitable activity in the same experimental conditions and energy range as 
the production data (see Fig, \25\-right panels) . Thus, very stringent acceptance 
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windows, which assure absence of any noise tail, can be considered and related 
efficiencies can be properly evaluated and used (see Fig. [26]) . In particular, to 
determine the efficiencies of the acceptance windows and their stability over 
the whole data taking periods, they are regularly measured by a suitable ^^^Am 
source ^^ I. These routine efficiency measurements are performed roughly each 
10 days; each time typically 10^ - 10^ events per keV are collected. 
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Fig. 26. Typical behavior of the overall efficiencies for the single-hit events. They 
are dominated by the efficiency related to the stringent acceptance windows used 
in the data analysis to assure the full rejection of noise near energy threshold (see 
Fig. [25]); see text. 



Let us now consider the radiation which can be induced in the crystal by 
the Dark Matter particles (see e.g. [^f5|6|8|T0l|TT] ) . Also thanks to the large 
number of available photoelectrons and to the single photoelectron threshold 
of each PMT working in coincidence, the detection efficiencies for internal 
photons and/or X-rays in the few keV range approaches the unity as well as 
for electrons. The detection efficiency for nuclear recoils is practically equal to 
unity considering the shortness of the recoil's range in matter. 

The energy threshold 

As shown, the new DAMA/LIBRA detectors have been calibrated down to 
the keV region (see Fig. [9] and Fig. [T9|) . This assures a clear knowledge of the 
"physical" energy threshold of the experiment. It obviously profits of the rela- 
tively high number of available photoelectrons/keV; the larger this number is, 
the lower is the reachable physical threshold. In particular, in DAMA/LIBRA 
the two PMTs of each detector work in coincidence with hardware threshold 
at single photoelectron level, while the software energy threshold used by the 
experiment is 2 keV. This latter value is well supported by the keV range cal- 
ibrations, by the features of the near-threshold-noise full rejection procedure 



^^ The use of ^^"^Cs Compton electron calibrations practically gives the same values 
as also already verified in DAMA/Nal. 
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and by the efficiencies when lowering the number of available photoelectrons 
(see for example Fig. [26]) . 



The single-hit scintillation events at low energy 

The procedure for noise rejection near energy threshold, described above, is 
the only procedure applied to the collected data. Fig. [27] shows, as example, 
the resulting cumulative low energy distribution of the single-hit scintillation 
events, as measured by DAMA/LIBRA detectors in an exposure of 0.53 ton x 
yr. This energy distribution is the mean value of all the used detectors; some 
differences among the detectors are present depending e.g. on their specific 
levels of residual contaminants and on their position in the matrix. 
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Fig. 27. Cumulative low-energy distribution of the single-hit scintillation events 
(that is each detector has all the others as veto) as measured by the used 
DAMA/LIBRA detectors in an exposure of about 0.53 ton x yr. The energy thresh- 
old of the experiment is at 2 keV and corrections for efficiencies are already applied. 
See text. 

It is worth noting that - as in the former DAMA/Nal experiment - neither 
other on-line or off-line techniques nor backgrounds from residual radioactiv- 
ity suppression are performed at all. In particular, background contribution to 
the counting rate of the single-hit scintillation events, arising from residual ra- 
dioactivity in the experimental apparatus, cannot precisely be extrapolated in 
the keV region e.g. because of: i) the limitation of MonteCarlo simulation pro- 
grams at very low energies; ii) the fact that often just upper limits for residual 
contaminants are available; iii) the unknown location of each residual con- 
taminant in each component; etc. Nevertheless, the investigations, presented 
in previous section, are extremely useful e.g. to qualify the detectors and to 
identify the sources which should be reduced in further developments of ra- 
diopure crystals, detectors' components, apparatus components, etc.. On the 
other hand, as known, the annual modulation signature, which is exploited by 
DAMA apparata, acts itself as an effective background rejection procedure. 

Response to nuclear recoils 

Finally, it is worth noting that, whenever WIMP (or WIMP-like) candidates 
are considered in corollary analyses of Dark Matter investigations, the response 
of the Nal(Tl) detectors to nuclear recoils has to be taken into account [4J. 
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Moreover, for a correct evaluation of the phenomenologies connected with the 
detection of WIMP (or WIMP-hke) candidates, the so-called Migdal effect [8] 
and the channeling effects [H] have also to be considered. Detailed discussions 
can be found e.g. in ref. 



Stability of the running conditions 

As already done for DAMA/Nal (see e.g. ref. [SE]), the stability of several 
parameters and of the running conditions is monitored along all the data 
taking. This aspect is directly related with the experimental results on the 
annual modulation signature and, therefore, it will be discussed elsewhere. 



9 Conclusions 



In this paper the ~ 250 kg highly radiopure Nal(Tl) DAMA/LIBRA appara- 
tus, running at the Gran Sasso National Laboratory of the I.N.F.N., has been 
described. The main components and performances have been discussed. 
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